Unusual linewidth dependence of coherent THz emission measured from intrinsic 
Josephson junction stacks in the hot-spot regime 
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We report on measurements of the linewidth A/ of THz radiation emitted from intrinsic Josephson 
junction stacks, using a Nb/AlN/NbN integrated receiver for detection. Previous resolution limited 
measurements indicated that A/ may be below 1 GHz - much smaller than expected from a purely 
cavity-induced synchronization. While at low bias we found A/ to be not smaller than ~ 500 MHz, 
at high bias, where a hotspot coexists with regions which are still superconducting, A/ turned out 
to be as narrow as 23 MHz. We attribute this to the hotspot acting as a synchronizing element. A/ 
decreases with increasing bath temperature, a behavior reminiscent of motional narrowing in NMR 
or ESR, but hard to explain in standard electrodynamic models of Josephson junctions. 

PACS numbers: 74.50.+r, 74.72.-h, 85.25.Cp 



Terahertz generation utilizing stacks of intrinsic 
Josephson junctions (IJJs) in the high-transition- 
temperaturc (high-T c ) cuprate Bi 2 Sr 2 CaCu208 + ,5 
(BSCCO) has been a hot topic in rece nt y ears, both in 
terms of experiment p]4Tl| and theory fl2H39l |. Typical 
emission frequencies f p are 0.4 - 1 THz, with a maximum 
output power of tens of //W emitted into free space ■ 
The IJJ stacks are usually patterned as mesas on top of 
BSCCO single crystals and contain 500 - 2000 IJJs. A 
large fraction of these IJJs oscillate coherently. There 
are good indications that the synchronization has to do 
with electric fields utilizing the mesa as a cavity: (i) f p 
increases inversely proportional to the junction width 
P, [l(| or fulfills the resonance condition for cylindrical 
mesas 0, E|, (ii) Low-temperature scanning- laser mi- 
croscopy imaged standing waves at bias points close to 
THz emission maxima 0, 0, H|. On the other hand, 
(iii) the linewidth A/ of THz radiation - so far mostly 
resolution limited to some GHz - is much lower than 
the estimated linewidth of the cavity resonances [jj HH 
and thus one cannot simply assume that the IJJs are 
slaved by the cavity; (iv) in some cases THz emission 
lines were fairly tunable or they occurred at positions 
not compatible with cavity resonances [IcJ. A recent 
mixing experiment [Toj has revealed a A/ as low as ~ 
0.5 GHz at low bias, rising even more questions about 
the mechanism of synchronization. 

In the IJJ stacks, operated at a bath temperature TJ, 
well below T c , there are two emission regimes. At mod- 
erate input power ( "low-bias regime" ) there is only little 
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heating, and the temperature distribution in the mesa 
is roughly homogeneous and close to Tt,. At high input 
power ( "high-bias regime" ) a hotspot (an area heated 
to above T c ) forms inside the mesa, leaving the "cold" 
part of the mesa for THz generation by the Josephson ef- 
fect. The hotspot seems to take part in synchronization 

To measure A/ and thus to gain more insight 
into the mechanism of synchronization, we have used 
a Nb/AlN/NbN superconducting integrated receiver 
(SIR)[40| for detecting THz emission from an IJJ mesa 
containing ~ 600 junctions, cf. Fig. [T] 

The 330 x 50 /im 2 wide and 0.9 /im thick mesa was 
patterned on a BSCCO single crystal (T c w 82 K), as de- 
scribed elsewhere @, 0, @ • Mesa and base crystal were 
contacted by Au wires fixed with silver paste. Contact 
resistances (~ 1 f2) have been subtracted in the data dis- 
cussed below. "Conventional" THz emission measure- 
ments were performed in Tsukuba, using a Fourier spec- 
trometer and a Si bolometer, as described in Mea- 
surements with the SIR were performed in Moscow. The 
SIR comprises on a single chip a planar antenna com- 
bined with a superconductor-insulator-superconductor 
(SIS) mixer, a superconducting flux-flow oscillator (FFO) 
acting as local oscillator (LO) and a SIS harmonic mixer 
(HM) for the FFO phase locking. The frequency of the 
LO is continuously tunable from 350 to 750 GHz, while 
the SIS mixer is effectively matched to the planar an- 
tenna between 450 and 700 GHz, limiting the SIR opera- 
tion range. A frequency resolution of the SIR well below 
100 kHz has been confirmed [4lj . 

Fig. [5] shows a typical set of data at Tf, = 40 K. The 
current- voltage curve (IVC) of the mesa is shown in Fig. 
[l^a). Upon increasing bias current J, near / ~ 25mA, 
one observes jumps (grey line of IVC) where the IJJs 
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Figure 1: (color online). Setup for detection of THz emission 
from intrinsic Josephson junctions (sketched at bottom). Top: 
central part of the superconducting integrated receiver (SIR) 
chip with a double-dipole antenna, a twin SIS-mixer and a 
harmonic mixer (HM) for flux-flow oscillator (FFO) phase- 
locking. For thermal shielding an infrared filter (2 mm thick 
GoreTex film) was installed between the BSCCO oscillator 
and the SIR. 



switch into their resistive state. When lowering I from 
~ 40 mA (black line) all junctions are resistive. Switch- 
back occurs below 3.65 mA. The outermost branch, where 
the emission data are taken, exhibits the typical back- 
bending arising from strong heating (33|. The high-bias 
regime embraces the region of negative differential resis- 
tance on this branch. Figs. [2£b),(c) show the broad- 
band THz emission power Pboh simultaneously detected 
with the IVC by a bolometer, vs. (b) current / and (c) 
voltage V across the mesa. In (b) one observes a broad 
(with respect to I) emission peak in the high-bias regime 
and a sharper one at low bias which, for this particular 
mesa, was more intense than the high-bias peak. Pboi 
vs. V, Fig. [He), shows a smooth peak for the high-bias 
regime, with a full-width-at-half-maximum in voltage of 
about 75 mV (black line) . By contrast the low-bias signal 
(grey line) seems to exhibit structure, indicating differ- 
ent groups of oscillating IJJs. Emission spectra, mea- 
sured using a Fourier spectrometer, are shown in Figs. 
H^d), (e). Spectrum (d) has been taken at point H in (a) 
where V = 736 mV, I = 19.43 mA. The emission peak is 
at f p = 618 GHz, corresponding to N = V/(/ p $ ) ~ 576 
oscillating IJJs. $o is the flux quantum. The (resolu- 
tion limited) lincwidth of the emission peak is ~ 15 GHz. 
Fig. [5] (e) shows a spectrum for bias point L where V 
— 697 mV, / = 3.76 mA. The emission peak, having a 
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Figure 2: (color online). Results of the THz emission experi- 
ment at 40 K: (a) IVC of the IJJ stack. THz emission power 
vs. current (b) and voltage (c). (d), (e) Fourier spectra of 
emitted radiation taken at points H and L in (a). Spectra of 
the mixing output of the SIR for (f ) three bias points in the 
high-bias regime and (g) for three bias points at low-bias. In 
(g) /lo was adjusted for each bias to bring the IF signal to 
the center of the spectrum analyzer, /lo values are indicated 
in the graph. Gray line indicates background noise. 



width which is also near the resolution limit, is at f p = 
594 GHz, corresponding to N ~ 567. 

The SIR measurements, showing the emitted radiation 
at the IF frequency, were performed near points H (Fig. 
H;f)) and L (Fig. ^g)). In Fig. [^f) the LO frequency 
of the SIR is /lo = 609.4 GHz. We show three measure- 
ments made for slightly (± 0.6 mV) different values of V. 
There is a single emission peak which shifts with increas- 
ing voltage. A/ of the two outermost peaks is ~ 60 MHz, 
while the center one is even sharper, A/ ~ 40 MHz. By 
contrast, the lincwidths seen in the low bias regime, cf. 
Fig. [U(g), ar c much larger and often exceeded the 6 GHz 
bandwidth of the IF amplifier. The lowest values are ~ 
0.5 GHz, in agreement with previous measurements [To| . 

Fig. [3] shows the most narrow line (A/ = 23 MHz) 
we have obtained in the high-bias regime. Data were 
taken at T b = 44 K, I = 16.3 mA and V = 0.721V. From 
/lo = 605.75 GHz and / p j F = 6.2 GHz we find f p = 
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Figure 3: (color online). Spectrum (most narrow emission line 
obtained) of the mixing output of the SIR for / = 16.3 mA, 
V = 0.721 V in the high-bias regime at T b = 44 K. The LO 
frequency is /lo = 605.75 GHz; / p ,if = 6.2 GHz. Dashed line 
is a Lorenzian fit with Af = 23 MHz. Inset shows drift of 
emission peak within 3 minutes. 



611.95 GHz and N 570. This value for Af is already 
quite close to the limit value that can be phase-locked by 
a regular room temperature semiconductor phase-locked- 
loop (PLL) system with regulation bandwidth of about 
20 MHz, presently limited by the delay in the cable be- 
tween the oscillator inside the cryostat and the PLL sys- 
tem (this limitation can be overridden by use of a cryo- 
genic PLL system with regulation bandwidth > 40 MHz 
[42] ) . We also investigated the stability of the THz emis- 
sion line. The inset of Fig. [3] shows two measurements, 
taken at a 3 min. interval. The 3 MHz drift observed was 
smooth and unidirectional and even lower than the typ- 
ical 3 min. drift (~ 10 MHz) of a free-running Nb based 
FFO. We also note that the emission lines observed at 
high bias have Lorentzian shape, cf. dashed line in Fig. 
[3j Recent 3D simulations [39[ based on the sine- Gordon 
equation and taking a hotspot into account yielded an 
asymmetric Fano shape in contrast to our experimental 
data. 

If Af of the high-bias THz emission signal were set by 
a cavity resonance, its quality factor Q = f j Af would 
have to amount some 10 4 . which is more than unrealistic 
for an IJJ stack 0, [3^, |43j]. At least for small mesas, 
cavity resonances become overdamped for T > 60 K [43j . 
For T ~ 40K we expect Q << 50. The low value of 
Af implies, that, when changing the bias current in Fig. 
[He), a sharp emission line "moves" through the AV = 
75 mV wide emission peak. The ratio V/ AV ~ 10 of this 
peak would indeed be compatible with the notion of a 
cavity resonance. 

In Fig. @] (a) we show by black circles Af vs. T},. 
Strikingly, Af decreases with increasing Xj, by about a 
factor 15 and, roughly, Af oc T b ~ 4 , cf. grey line in Fig. 
Eta). For comparison we have also added for two temper- 
atures low-bias data for f p near 616 GHz. In this regime 
it was difficult to obtain systematic data points, because 
at most temperatures only a broadband (> 6 GHz) signal 
appeared in the SIR output. This further indicates that 
at low bias phase-lock is incomplete. Noting that Tj, may 
not be the best variable for our system, in Fig. 0] (b) we 
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Figure 4: (color online). Linewidth vs. (a) bath temperature 
and (b) dc input power. In (a) solid (open) circles are for high 
(low) bias. In (b) open squares are for low bias, other symbols 
for high bias. Solid (black) circles: same data as in (a). Solid 
(blue) squares: linewidth at T;, = 40±1K and various values 
of f p . Linewidth of Fig. [3]is shown by solid (green) diamond. 



also show Af vs. the dc input power Pj c = IV, where 
for V we have used the measured value, not corrected 
for contact resistance. Low-bias data are shown by open 
squares. The black circles show the same high-bias data 
as Fig. Eta); Af increases with increasing Pdc- Note that 
the smallest Pd c is at high Tb and the highest Pj c is at the 
lowest Tb- This does, however, not imply that the sam- 
ple is hotter at low Tb- Simulations of the temperature 
distribution, performed analogous to the calculations in 
[33|, yielded temperatures of the cold part of the mesa 
which arc some 10-15 K above but still increase with 
increasing Tj,. For the hot part we find peak temper- 
atures Th in the range 80-110 K, with the highest value 
indeed at the lowest value of Tb. For comparison we show 
by solid squares in Fig. |4jb) Af out of a measurement 
series, where we have fixed Tb to 40±1K and varied f p . 
Within this curve f p varies monotonically from 697 GHz 
to 562 GHz. In contrast to the f p = const, curve, Af 
is about constant at intermediate values of Pj c and in- 
creases outside of this regime. Thus there seems to be no 
systematic dependence of Af on Pj c - Further note, that 
a slight increase of Af at low Pd c is visible both on the 
Tb = const, and the fp = const, curve. Here, the bias 
is close to the hotspot nucleation point and it is possible 
that the hot part of the stack already is close to or even 
below T c , i.e. we approach the low-bias regime. 

We next compare the high-bias data to theoretical pre- 
dictions. We decompose Af into a sum of two contribu- 
tions A fx and Af x , as it has been done in an analysis 
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of the radiation linewidth of the FFO (44[. A/t arises 
from fluctuations of the quasiparticle current (Nyquist 
noise) and Af x represents additional fluctuations of - 
at this point - unspecified origin. For a current biased 
lumped junction Af T = (47r / ! §l){R 2 d / 'R)k B T ', where R is 
the junction resistance, Rd is the differential resistance of 
the IVC at the bias point and T is the junction tempera- 
ture [mm)]. For a phase-locked ID array of N junctions 
the linewidth is expected to decrease oc TV -1 (46[. In 
our case, in the high-bias regime the hot part of the IJJ 
stack is an excellent candidate for a shunting network 
causing phase- locking [1, [j| HH, H|| . Our IJJ stacks are 
not lumped and also T not only differs from T}, but also 
varies strongly within the stack. Nonetheless an anal- 
ysis of Afx may give some indication whether or not 
A/ vs. T\y can be explained via Nyquist noise. We first 
look at the cold part of the mesa only. The out-of-planc 
resistance R c of BSCCO has a negative temperature co- 
efficient (roughly R c oc T^ 1 ) and thus T/R oc T 2 . We 
cannot infer Rd from the IVC, since T varies strongly 
along this curve. However, from small mesas it is known 
that Rd is smaller but not very different from R. Thus, 
an unrealistic R d oc T~ 3 would be required to explain 
the data. We also considered a hot resistor (resistance 
R s , temperature T s ), representing the hotspot, in parallel 
to a cold one (resistance Rj, temperature Tj), as it has 
been analyzed in [47j for a single junction, with the re- 
sult that R should be replaced by R e g = (Rj 1 + R^ 1 )^ 1 
and T by T cff = T s /R s + Tj/Rj. R cS follows from the 
measured IVC. With T s ^T h » Tj and R s « Rj we 
get T e ff ss Th and, with Rd ~ i? e ff, one obtains an almost 
temperature independent Afx ~ 250 MHz. Dividing this 
by N yields a linewidth which not only has the wrong 
temperature dependence but also is much lower than the 
measured A/. 

The second contribution to A/, Af x , may arise from 
moving fluxons. This is the case for the FFO where, 
however, Af x oc T 1 / 2 was found 0. Recently, THz 
generation in IJJ stacks by bound fluxon-antifluxon pairs 
(breathers) oscillating near a cavity resonance was pro- 
posed. Whether or not this mechanism leads to the 
observed T depencence of A/ is unclear, although we 
would expect an increase with T, as for the case of the 
FFO. Also, more conventional simulations based on 1D- 
coupled-sine-Gordon equations yielded a Af which in- 
creased with temperature [H| . 

Another hint may be, that at low Tj, the in-plane ther- 
mal gradient in the stack is larger than at high T\,. How- 
ever, we cannot translate this to the observed Af vs. Tf, 



dependence. 

Let us thus also look at other effects in physics that 
cause a decrease of Af with increasing temperature. Mo- 
tional narrowing in NMR or ESR are prominent examples 
[H, 51] • R is tempting to speculate about an analogous 
effect in the IJJ stacks. In a nutshell, motional narrow- 
ing occurs when some spatially varying field, causing line 
broadening, is reduced on average due to temporal fluctu- 
ations. In our case this could be due to a time-dependent 
hotspot. Generally, it is possible that the hotspot oscil- 
lates We have searched for such oscillations for fre- 
quencies in the kHz and MHz range but could not find 
any. Still, the temperatures of the hot and cold parts 
of the stack may fluctuate, causing fluctuations of the 
hotspot edge position and as a result also variations in 
junction parameters such as critical current and resis- 
tance. Suppose that the static spread in these parameters 
is too large to achieve complete phase-lock. Temporally 
fluctuating parameters could lead to a situation, where 
for some time these parameters are in the synchroniza- 
tion window, causing improved phase-locking. 

To conclude: At low bias the best values for Af are 
on the order of 0.5 GHz, consistent with previous mea- 
surements pfl| . However, the irregular lineshape and 
the strongly varying linewidth indicate that not all junc- 
tions are phase-locked and a superposition of several lines 
is observed. In the high-bias regime we observe much 
lower values of Af down to ~23MHz. The improved 
linewidth compared to the low-bias regime suggests that 
the hotspot acts as a shunting element providing phase- 
lock. Af, measured at fixed emission frequency, de- 
creases with increasing temperature, reminiscent on mo- 
tional narrowing observed in NMR or ESR. This tem- 
perature dependence is hard to explain within standard 
models for the linewidth of oscillating Josephson junc- 
tions. In any case, the narrow emission line observed in 
the high-bias regime makes intrinsic Josephson junction 
stacks interesting for, e.g., THz spectroscopy. 
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